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Abstract. The revision of the photospheric abundances proferred by As- 
plund et al has rendered opacity theory inconsistent with the seismologically 
determined opacity through the Sun. This highlights the need for a direct seis- 
mological measurement of solar abundances. Here we describe the technique 
used to measure abundances with seismology, examine our ability to detect dif- 
ferences between solar models using this technique, and discuss its application 
in the Sun. 



Introduction 



A decade ago, solar models matched the Sun remarkably well. Sound speed 
profiles in models differed from the Sun's actual sound speed by less than half a 
percent. Then Asplund et al. | ( |2005 ) reanalyzed the solar optical spectrum using 
improved atomic physics and a 3D hydrodynamical model of the atmosphere 
(instead of a ID hydrostatic model). Their results indicate that the abundances 
of the heavy elements should be lower by a substantial amount. For example, 
carbon was lowered by 35%, nitrogen by 27.5%, oxygen by 48%, and neon by 
74%. Solar models that are evolved using the new abundance mixture give worse 
agreement with helioseismic constraints than models using the old abundances. 
Problems include a tripling in the sound-speed discrepency below the convection 
zone, a convection zone that is too shallow, and a helium abundance in the 
convection zone that is too low. 

Many attempts have been made to reconcile the new abundances with seis- 
mology. A few examples include increasing the opacities below the convec- 
tion zone, increasing abundances within their uncertainties, increasing the neon 
abundance, enhancing diffusive settling rates, and including early accretion of 
lower-Z material. Various combinations of these changes have also been explored. 
None of the alterations has resolved the discrepency between the revised abun- 
dances and se ismology. Fo r furt her discussion of the problem and the attempted 
solutions, see Basu (20081 and Guzik et al. (20081 in these proceedings. Since 



the mismatch between seismology and the new abundances remains a problem, 
a seismic determination of the heavy element abundances in the solar convection 
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zone is necessary. The techniques for a seismic measurement of heavy elements 
are discussed in the following sections. 



2. Techniques 



Dappen & Gough (1984) measured the solar helium abundance through the 



effect of helium ionization on the adiabatic exponent 71. Ionization lowers 71 
from the ideal gas value of 5/3. Since the convection zone stratification is almost 
purely adiabatic, the gradient of the square of the sound speed in that region 
can be written as 
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This expression can be separated into a term called W that contains the seismic 
variables and a term called that contains only thermodynamic variables. 
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Since contains a derivative of 71, it shows a more pronounced variation due to 
the ionization of each element than 71 alone shows. For the helium abundance 
measurement, Dappen and Gough used seismic data to obtain the sound speed, 
then computed W and for their model. Figure [T] shows the helium hump that 
Dappen and Gough discovered in W and where He II ionization takes place 
in the convection zone. Here W and are ploted against the acoustic radius 
t = J c dr. The signatures of heavy-element ionization can already be seen in 
the small wiggles lower in the convection zone (Gough 2006). Near r = 0.63, 
we see C V ionization. From 0.45-0.55, O VII ionization appears. Near 0.58, 
N VI shows up as an imperceptibly small hump that is impossible to measure 
here. Near 0.68, there is a double hump due to the first two K-shell ionizations 
of Ne. The larger humps are merged ionization zones. Since the humps from 
heavy element ionization are such small features, great care must be taken to 
obtain an accurate measurement of the abundances. 



2.1. Inversion 

It is common to use sound speed and density or sound speed and 71 as solar 
inversion variables. However, the extremely small ionization humps of the heavy 
elements require a more sensitive inversion. We suggest inverting the frequency 
differences directly for W and 0. 
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(K W>Q SW + K etW 50) dr 
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In the adiabatically stratified convection zone, W = 0. This allows us to rewrite 
the integral as 

5to 2 r R 

-2 = / (K WtQ + K e ,w) SQ dr. (5) 
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Figure 1. 



W and 6 throughout the convection zone in model S of 
(19961. The large hump around r = 0.82 is 



Christensen-Dalsgaard et al 



due to helium ionization. The smaller wiggles lower in the convection zone 
are due to heavy clement ionization. 



Since this formulation has the added advantage of reducing the problem to a 
single-variable inversion, we no longer need to worry about contamination from 
a second inversion variable. 

2.2. Isothermal Sound Speed 

The conversion from the familiar sound speed and density kernels (see, for ex- 
ample, Gough Thompson |199l" ) to the desired W and G kernels requires quite 



a bit of algebra. One possible trick to simplify the calculation is to reconfigure 
W and in terms of the square of the isothermal sound speed u = p/p. The 
gradient of u in an adiabatically stratified region is 



du Gm ( 1 



dr r 2 \7i 



(6) 



This leads to a set of definitions analogous to equations ^ and ^ for the 
separation of seismic and thermodynamic variables: 

Gm dr 
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The derivative of 71 in the definition of does not appear in the new 0. This 
greatly simplifies the kernel conversion. However, without the 71 derivative, O is 
less sensitive to the effects of ionization than is. As shown in figure [2J O does 
not have visible heavy-element ionization humps. Since the effect of ionization 
is so small to begin with, reducing it even further by switching to is not worth 
the gain in algebraic simplicity in the kernel conversion. The 71 derivative is 
essential to highlight the effect of heavy-element ionization. We therefore choose 
to use the original definitions of W and G. 
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Figure 2. 6 and throughout the convection zone in model S of 



Christensen-Dalsgaard et al. (19961. The heavy-element ionization humps 



are clearly visible in 0, yet they do not appear in 0. 



2.3. Pressure Ionization 

One source of uncertainty in solar equations of state is the effect of pressure ion- 
ization. Pressure ionization occurs when matter is so dense that electron orbits 
overlap. The electrons in the overlapping orbits no longer belong to a partic- 
ular atomic nucleus. The pressure effectively frees these electrons, increasing 



the ionization of the plasma. Baturin et al. (20001 developed a technique to 



account for pressure ionization in an equation of state by applying constraints 
to the effective atomic and ionic sizes. They used a single pressure-ionization 
parameter £ for all species of ions. They then implemented a model in which 
they calibrated ( to fit solar data. The calibrated pressure-ionization model 
fit the data better than the parameter-free equations of state, highlighting the 
importance of pressure ionization for an accurate equation of state. 
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Changing the effective atomic size in the pressure-ionization model leads 
to a vertical shift in 71. Since 71 appears as a factor in 0, 6 also suffers from 
this shift due to the atomic size parameterization. Therefore, we cannot just 
measure the average value of to obtain the heavy-element abundances. In 
order to obtain an accurate seismological measurement, we must consider the 
effective atomic size by using a pressure-ionization model for the equation of 
state or measure the humps relative to the background level of 6. 

3. Conclusions 

In order to obtain a seismic measurement of the heavy-element abundances in 
the Sun the following techniques should be used: the frequency data should be 
inverted directly for W or instead of sound speed, the original definitions for W 
and O should be used instead of the W and derived from the isothermal sound 
speed formulation, and the effect of atomic size calibration on the background 
level of should be taken into account. These choices are crucial for accurately 
measuring the small effect of heavy-element ionization in order to determine 
solar abundances seismically. 
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